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Abstract 

I analyze the Nonlocal Conserved Kardar-Parisi-Zhang (NCKPZ) equation with spatially cor- 
related noise. This equation is also known as the Nonlocal Molecular Beam Epitaxy (NMBE) 
equation andv was originally suggested to study the effect of the long-range nature of interactions 
coupled with spatially correlated noise on the dynamics of a volume conserving surface. I find an 
exact result for a subfamily of NCKPZ models in one dimension, and thus establish an exact result 
for MBE processes for the first time. Then, to complete the picture, I construct a Self-Consistent 
Expansion (SCE) and get results that are consistent with the exact result in one dimension. I 
conclude by discussing previous results obtained for NCKPZ using dynamic renormalization group 
approach, and find that this approach does not yield the exact result I obtain. This discussion 
shows the advantage of the using SCE to deal with non-linear stochastic equations. 
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The field of disorderly surface growth has received much attention during the last two 
decades. A multitude of different phenomena such as fluid flow in porous media, propagation 
of flame fronts, flux lines in superconductors not to mention deposition processes, bacterial 
growth and "DNA walk" are all said to be related to the famous Kardar-Parisi- Zhang 
(KPZ) equation 

In spite of the great success of KPZ to describe many growth models and phenomena, 
there has been a growing pool of data that is not well described by KPZ, that called for 
further research. One of the first classes that belongs to this non KPZ behavior is the well 
known conserved KPZ equation class (sometimes also called the Molecular-Beam-Epitaxy 
(MBE) class) j^l-j^l- However, the modified behavior introduced by the MBE equation is 
not at all sufficient to account for all the rich non KPZ experimental data in the field. 

A different line of research suggested that the basic growth equation should not be 
changed. Instead, the white noise that app ears in the original KPZ equation should be 
either temporally or spatially correlated |lfll|-jl^. This approach was indeed quite success- 
ful, but it still failed to give a good account for all the measured scaling exponents. 

Recently, Mukherji and Bhattacharjee [l^J proposed a phenomenological equation that 
takes into account the presence of long-range interactions. This equation, known as the 
Nonlocal KPZ (NKPZ) equation, has a nonlinear term that couples the gradients of the 
growing surface at two different points. This equation was studied using the Dynamic 
Renormalization Group (DRG) method [13] and so was its generalization that includes 
spatially correlated noise Q]. It turned out that a very rich family of phases was found. 

Just like in the local case, a nonlocal conserved KPZ (NCKPZ) equation was needed in 
order to account for the dynamics of a volume conserving surface that experiences long-range 



interactions. This equation was proposed by Jung et al. [15|,^^. The continuum equation 
they suggest for the height of the surface h (r, t) at a point r and time t measured relative 
to its spatial average is 

dh (r, t) 



dt 

xVh{f+f',t)-Vh{f-f^,t) + 7]{f,t), (1) 

where K is a constant, the kernel g (r) represents the long-range interactions with a short- 
range part Aq^'' (r) and a long-range part ~ Xpr''~'^, where d is the substrate dimension and 
p is an exponent characterizing the decay of the long range interaction. More precisely in 



Fourier space g (g) = Aq + XpQ ^. Finally, r] (r, t) is a conservative spatially correlated noise 
term that satisfies 



{Vir,t)) = 

{r] (f, t) r] (r', t)) ~ |f - r']" 5 (t - t') , (2) 

where a is an exponent characterizing the decay of spatial correlations (in the case of con- 
servative white noise , which is of great interest in the MBE system, a = 0). 

After suggesting the model, Jung et al. ll^ investigated the NCKPZ problem using 
DRG analysis, and derived a very rich and complex picture of phases - each one of them 
described by a different set of critical exponents. In this letter I will not go into the details 
of all those phases, but would like to focus on the strong coupling behavior of NCKPZ. In 
this regime the DRG analysis yields the power counting result 

a= (2-d-p + 2a)/3, z = (rf + 10 - 2p - 2a)/3, (3) 

where a is the roughness exponent and z is the dynamic exponent. This solution is said to 
valid for positive values of p and ford<2 + 2p + 2cr. 

In this letter I will show that the above DRG analysis does not give the full description of 
the different phases in the strong coupling regime of NCKPZ. More specifically, I will show 
that an exact solution is possible for some one-dimensional cases (namely for p = 2 + 2a). It 
is interesting to mention that this is the first time an exact result is suggested for the family 
of MBE processes. It will then be immediately seen that this exact solution is not compatible 
with the DRG tesuh. Thenjn orietto citify the situation, I will apply a method developed 



by Schwartz and Edwards 



la 



22l l23l | (also known as the Self-Consistent-Expansion (SCE) 
approach). This method has been previously applied successfully to the KPZ equation. The 
method gained much credit by being able to give a sensible prediction for the KPZ critical 
exponents in the strong coupling phase, where, as is well known, many Renormalizat ion- 
Group (RG) approaches failed (as well as DRG of course) j^]. For the specific problem of 
NCKPZ, the SCE method indicates that two strong coupling solutions exist, one of which 
agrees with the power counting solution also obtained by DRG, while the other one recovers 
the exact one-dimensional result (and also extends it to higher dimensions). It turns out that 
each solution is valid for different values of the parameters cr, p and d (mutually excluding 
values) . 



In order to derive the exact result, I have to consider eqs. (0)-(|21) in Fourier components 



dh 1 

-Kq^hq ^ (^0 + Apg^^j (i- m) Sq^i+rnhih^ + rjq, (4) 



where rjq is the Fourier component of the noise that satisfies: 

= 

it) W it')) = 2Doq'-'-5q^q,6 {t - t') . (5) 

I will show that for the case Aq = an exact solution is possible for the steady state 
distribution resulting from eq. Q. However, before doing that I would like to mention that 
when p > 0, the Ap term in eq. dominates over the Aq term in the asymptotic small 
q limit, so that this term can be neglected. On the other hand when p < the Ap term 
can be neglected compared to the Aq term, and I retain the standard MBE results, unless 
Aq = exactly (The specific case Aq = and p < is reasonable only under special physical 



conditions j25|). It will be understood therefore, that while the results obtained for positive 
p's hold also for Aq 7^ 0, the results obtained for negative p's hold only when Aq is strictly 
zero. In other cases (i.e. p < and Aq 7^ 0) the MBE results are recovered. 

I now turn from the NCKPZ equation in Langevin form to a Fokker-Planck form 



dP 
'dt 



^dhq 



Doq 



2~2a 



d 



dh. 



+ Kq%+ 



A.a^-P 



(i- mj 6q^i+rahtK 



P = 0, 



(6) 



(where P ({/ig}, t) is the probability of having the configuration {/ig} at a specific time t). 

Consider now the family of models defined by p = 2 + 2a. By direct substitution, it can 
be verified immediately that the following probability function 



Fo = A^exp 



h h 

q r'q 



(7) 



(where N is the normalization yielding / Po{h} Dh = 1, and (pq is the two point correlation 
function) is a solution of the steady-state Fokker-Plack equation (i.e. with dP /dt = 0) for 
the one-dimensional case [d = 1) if (pq = Aq~^ with A = Do/ K and F = p. Translating this 
result to the frequently used notation in the field of surface growth I get 



" = (8) 

where a is just the roughness exponent. This result is therefore consistent with the known 
exact solution of one-dimensional KPZ a = 1/2 (by setting a = and p = 2). However, this 
result is manifestly different from the DRG result cited in eq. (jH))- 

After establishing the exact result I turn to the application of the SCE method to NCKPZ 
in order to clarify the conflict as well as to complete the picture for any noise exponent 
a and for higher dimensions. The SCE method is based on the Fokker-Planck form of 
NCKPZ equation, which is the reasonable thing to do once you have an exact solution for 
that equation in one dimension. The method constructs a self-consistent expansion of the 
distribution of the field concerned. The expansion is formulated in terms of (pq and Ug, where 
(pg is the two-point function in momentum space, defined by (/)g = {hgh_g)g, (the subscript 
S denotes steady state averaging), and Ug is the characteristic frequency associated with hg. 
I expect that for small enough q, (pg and Ug obey power laws in q, 

(f)g = Aq~'^ and Ug = Bq^ . (9) 

[Since dynamic surface growth is a remarkably multidisciplinary field, there are almost 
as many notations as there are workers in the field. Therefore I give a brief translation of 
our notations to those most frequently used: 

z = iji^ a = {T-d)l2 and /? = (P - rf)/2/x.] (10) 

The SCE method produces, to second order in this expansion, two nonlinear coupled 
integral equations in (pg and uog, that can be solved exactly in the asymptotic small q limit 
to yield the required scaling exponents governing the steady state behavior and the time 
evolution. 

By defining Mgi^ = ^'f"'' if' ^q,i+m, Kg = Kq'^ and Dg = Dnq^'^ it can be seen 
that the NCKPZ equation is of the general form discussed in refs. 
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where the SCE 



method is derived. Working to second order in the expansion, this yields the two coupled 
non-linear integral equations 



Dg - KgCPg + 2j2 



Mqim Mgim (pe(pm 



5 



_ 2 ^ ^qim Mirnq 4>m4'q _ 2 ^ ^ qlm^ miq4> l4> q 



£,m 



UJ„ + Ui + UJr, 



and 



UJn 



Mimq4>m + Mmlq4>i 



qim ' 



(11) 



(12) 



where in deriving the last equation I have used the Herring consistency equation 26] . In fact 

eading to a different consistency 
23|, that this does not affect the 



Herring's definition of is one of many possibihties, each 
equation. But it can be shown, as previously done in ref. 
exponents (universality). 

A full solution of equations (fTTjl and (fT^ in the limit of small q's (i.e. large scales) yields 
a very rich family of solutions that will not interest us here, but may be the subject of a 
future detailed paper. Instead, I focus on the strong coupling solutions obtained by SCE. As 
mentioned above, two strong coupling solutions are obtained by SCE, each of them describes 
the critical exponents for a different set of the parameters p, a and d. 

The first solution is obtained just by power counting, and reads /i = (c? + 10 — 2p — 2cr)/3 
and r = (rf + 4 — 2p + 4cr)/3. This solution is valid for p — 1 < a, (i < 6 — 2cr + p, 
(i < 2 + 2cr + 2p and < 8 + 2cr — 4p. Therefore, by direct inspection, this solution is 
not valid for the one-dimensional case with p = 2 + 2(T, and hence does not recover the 
exact Gaussian result obtained above. It should also be mentioned that when translated to 
frequently used notation (using eq. (fTUI) ) it can be seen that this solution is actually the 
same as the one obtained by the DRG analysis as in eq. 0- 

The second solution is a non power counting solution. It turns out that for d = 1 eq. 
flllll is exactly solvable, and yields T = p. This result is actually the exact Gaussian solution 
of the one- dimensional case, mentioned above. In addition, from eq. (fT^ . the dynamical 
exponent can be extracted for a certain range of the parameters, namely 1 < p < 3 and 
p > (4(7 + 5), where it becomes = (9 — 3p). 

For dimensionality higher than one (i.e. d > 2) such an exact solution in closed form 
cannot be found, and the strong coupling solution is determined from the combination of 
the scaling relation n = {d + 8 — T — 2p)/2, and the transcendental equation F (F, p) = 0, 
where F is given by 



F (F, p) = - / dH 



t- (e - 1) 



e-t 







[{e.t) 


e-i 



' r ^ + e ■ (e - t) t 



2-p 



e-t 



+ 



6 




ti'+ e-t + 1 



t- (e - tj 



t 



-r 



e-t 



(13) 



and e is a unit vector in an arbitrary direction. The last equation has to be solved numerically 
of course. It addition, This solution is valid as long as the solutions of the last equations 
satisfy the following four conditions: (i + 4 — SF — 2p + 4cr < 0, d — 3T + 2p < 0, d <T and 
d-A-r + 2p <0. 

It should be mentioned that SCE yields many other "weak coupling" solutions to NCKPZ. 
They are called "weak coupling" in the sense that all of them yield the exponents of the 
linear MBE problem with correlated noise in the limit of p = 0. These solutions are also 
not obtained by DRG. 

In summary, I have established an exact solution for the one-dimensional NCKPZ equa- 
tion in the special case where p = 2 + 2a. Then, I calculated the critical exponents for the 
most general case (i.e. any p, a and d) using the SCE method, and obtained two possible 
strong coupling solutions. The first strong coupling solution is a power counting one, and is 
identical to the previous result of the DRG analysis. The second strong coupling solution 
is non power counting, and is not obtained by DRG. For d = 1, this solution turns out to 
be the exact Gaussian solution obtained at the beginning. This suggest that there is a big 
advantage in using the SCE to deal with such nonlinear problems. 

Another advantage of using SCE has to do with the well known technical problem raised 
by Janssen j2l|. Janssen suggested that there should be a correction to the scaling relation 
because of a nontrivial renormalization of the coupling constant. This correction appears 
in a two-loop RG calculation, but not in a one-loop calculation. Janssen himself suspected 
that in order to avoid this difficulty, a mode-coupling approach might be more appropriate. 
The SCE method is indeed closely related to the mode-coupling approaches, in the sense 
that similar (but not identical) equations are obtained, while the underlying derivation is 
different. Therefore, the fact SCE recovers the one-loop DRG power counting results of Jung 
et al. , for a specific range of the parameters Q| , is a nontrivial corroboration of that 
result. 
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